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^H ■ ABSTRACT 

(^: 

^i^ , Clumped dark matter arises naturally within the framwork of generic cosmological dark 

tG ' matter models. Invoking the existence of dark matter clumps can also solve may unexplained 

^~^, mysteries in astrophysics and geology or geophysics , eg. the galactic gamma-ray halo and 

Q ■ the periodic terrestrial flood basalt volcanic episodes. Clumped dark matter is dynamically 

M I stable to friction and will not heat the disk. Such clumps may have already been discovered 

^ ■ in the form of dwarf spheroidals, and further searches are encouraged by the results of this 

C^ ■ 

. . , paper. 

>. 
rS I 1 Clumping as a Prediction of Cosmological Models 

P.: 

As the bulk of the matter in the Universe is dark, a better understanding through 
theoretical models is important, and much current activity is taking place in that 
direction. Undoubtedly more significant than this would be the detection of this dark 
matter. Much progress is taking place, and it is possible that dark matter may have 
been already detected. A Gran Sasso group claims to have found anomalous signs 
of a 56 GeV dark matter candidate fll . For recent discussions of the DAMA results, 
consult [0, PI- However, a definite confirmation of this is awaited. 

As the visible matter clumps together to form stars, planets, etc. an interesting 
question is whether the dark matter also displays this tendency of clumping. In- 
terestingly several dark matter models do suggest that clumps of dark matter arise 
naturally during the course of evolution of the universe. Silk and Stebbins [^ consid- 
ered cold dark matter models with cosmic strings and textures appropriate for galaxy 
formation. They found that a fraction 10"'^ of the galactic halo dark matter may exist 
in the form of dense cores. These may survive up to mass scales of 10®Mq in galaxy 



halos and globular clusters [^]. Analysing the stability of these clumps of dark mat- 
ter, they found that the cores of these clumps will not be affected, although the outer 
layers may be stripped off by tidal forces. In the cosmic string model, the dumpiness 
C, defined as the ratio of clumped matter concentration to normal concentration, of 
dark matter at the present epoch would be @] 

C ~ IQ'^fcih^^l (1) 

where fd is the fraction of dark matter in clumps, H is the Hubble parameter 
parametrised as 100 h km/ s Mpc~^, and Qq is the closure energy density of the 
Universe. 

Subsequently Kolb and Tkachev studied isothermal fluctuations in the dark 
matter density during the early universe. If the density of the isothermal dark matter 
fluctuation or clumps, $ = Spdm/pdm, exceeds unity, a fluctuation collapses in the 
radiation-dominated epoch and produces a dense dark matter object. They found 
the final density of the virialized object pp to be 

Pf ~ 140<l>3(<l> + l)peg (2) 

where peq is equilibrium density. 

For axions, a putative dark matter particle, density fluctuations can be very high, 
possibly spanning the range 1 < $ < 10^. The resultant density in miniclusters can 
be as much as 10^° times larger than the local galactic halo density. The probablility 
at present of an encounter of the earth with such an axion minicluster is 1 per 10^ 
years with $ = 1. Kolb and Tkachev found two types of axion clumps arising from 
two kinds of initial perturbations: 

• Fluctuations with 10~^ < $ < 1 collapse in the matter-dominated epoch. 

• Fluctuations with $ > 1 collapse in the radiation-dominated epoch. 

If the dark halo is mostly made of neutralinos, then the clumping factor in the 
MSSM could be less than 10^ for all neutralino masses p. 

Structure models involving the hierarchical accretion and mergeer of dark mat- 
ter halos are well-motivated and attractive. Numerical simulations of the substruc- 
ture with galactic and cluster halos which form within hierarchical universes are 
easily reproduced, yielding thousands of substructure clumps. The models analsed 
by B.Moore, S.Ghighna et al. [0 predict that the Milky Way should contain at least 
500 satellites of bound masses > IO^Mq and tidally limited sizes > kpc. Remarkably, 
dark matter substructure survives on galactic scales, so that galaxy halos appear as 
scaled versions of galaxy clusters. One model considered by S.Ghigna, Ben Moore 
et al. 1^ resolved 150 halos within a model galaxy cluster. Sazhin et al. [^ have 
suggested that microlensing events observed towards the Large Magellanic Cloud may 
consist of dark matter clumps. 



2 Astrophysical Constraints and Stability of Clumps 

Astrophysical constraints can be used to place bounds on the mass of the clumps. 
The clumps are also stable with regard to main astrophysical disruption processes. 

2.1 Disk 'Puffing' by Dark Clumps 

As the dark matter clumps traverse the Galactic halo, they will impart energy to the 
stars residing there. This leads to a gradual puffing up of the disk, and the stars 
present are heated. In this process, the older stars get heated more than the younger 
ones. Lacy and Ostriker considered the case of black holes of mass lO^M© generating 
the observed puffing of the galaxy, and concluded that this was the best mechanism to 
explain the observed amount of 'puffing' of the galactic disk |l^. However, subsequent 



data revealed that the velocity dispersion of disk stars, a, may no longer rise as quickly 
as ^/t. which was one of Lacy and Ostricker's assumptions. In addition, other sources 
of heat like spiral density waves and giant molecular clouds may give a better fit to the 
data. Hence the Lacy-Ostriker model of black holes generating the observed puffing 
of the disk is not as convincing now as it was when first suggested. However, they 
obtained an upper limit on the density of halo objects of mass M, which still hold : 

Qb < ^h min [l, {M/Mheat)-'] (3) 

where, hheat = 3 x 10^ (i^/lO^^)"' (4) 

where tg is the age of the galaxy. This condition must hold, otherwise the disk would 
be more puffed up than observed. Using this formula, one can obtain an upper limit 
on the mass of the halo objects in any galaxy. Applying it to the Milky Way yields an 
upper limit on the mass of dark matter clumps of M < 2 x lO^M©. Silk and Stebbins 
obtained, using slightly different arguments, a similar limit of lO^M© to avoid the 
problem of unacceptably heating the disk Q, the general limit for halo objects so as 
not to heat the disk. However, if the halo objects form pregalactically, which is the 
case if they are the precusors to the galaxy, then halo objects should have the same 
order of mass. Thus from an analysis of the gas-rich dwarf galaxy DD0154 one finds 
M < 7 X IO^Mq. For the dwarf galaxy GR8 the limit is M < 6 x lO^M©. 

Disk heating can also lead to 'gravitational shocks'. During orbital passage of a 
dark matter clump through the inner galaxy, the clump is heated. This is a major 
factor in the disruption of globular star cluster and hence is expected to be a primary 
factor in the disruption of cold dark matter clumps |ll|]. This process is considered 
in the next section. 

2.2 Dynamical Friction 

Gravitational shocking leads to considerable destruction of clumps. The destruction 
timescale for clumps is given by |TT| 



tdest ^ 4 X ;Lg4 / lOpC \ f Rperi \ ( Vclump \ ( 10^°Mq \ ( M^ump X /^x 



\{r clump)} \Skpc) \imkm/s) \Mi^ige{Rperi)l VIO'^0/ 



where Mhuige is the spherical or disk mass interior to the clump orbit, 

< fdump > is the mean clump radius, 

t orbit is the internal orbit time,, 

tdest is the destruction time scale, 

Mq is the solar mass, 

Mciump is the mass of a clump, 

Vclump is the mean clump velocity, 

Rperi is the perigalactic distanve of the clump orbit. 



Numerical simulations give a clump destruction rate of ||Tl| 

4:r - f #^) '^' (-^^] ''' ^Istr (6) 

\ -^''^ clump / \ Pclump / 

where i/^g^^ is the present globular cluster destruction rate, pgc is the density of 
glonular clusters and Mgc their mass. 

The present globular destruction rate is obtained from numerical simulations in- 
cluding the effect of gravitational shocking by the disk and bulge and evaporation 
losses. From such calculations it is inferred that this rate is 10~^^y~^. 

Thus, if the mass of the clump rriciump > 10'^ Mq, clump destruction is unimportant 
via dynamical processes over the last 10^° years 0]. Clumps in the inner few parsec 
of the halo are likely to have been destroyed, while clump halos in the outer part of 
the galaxy are likely to have survived intact. If a fraction of / ~ 1% of the halo has 
survived in clumps, and if Pdump ~ ^^^Phaio, then the resultant 7 ray signal leads to 
10 times larger fluxes than for a uniform dark matter halo. 

Tormen, Diaferio and Syer |T^ also analysed the survival of substructure in dark 
matter halos. Using high-resolution N-body cosmological simulations, they analysed 
the survival of dark matter satellites falling into larger bodies. They found that all 
satellites preserve their identity for some time after merger. Satellites with less than 
a fre percent mass of the halo may survive for several billion years, whereas larger 
satellites rapidly sink into the center of the main halo potential well and soon lose 
their identity. 

2.3 Diffusion Towards the Galactic Center 

As halo objects traverse through the halo, they will lose energy and consequently 
drift towards the Galactic nucleus. It has been shown |T0| that halo objects will be 
dragged into the nucleus by the dynamical friction of the spherical stars from within 
a galactic radius 
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Rdf = (m/IO^Mq) (t,/10^^) kpc (7) 

and that the total mass dragged into the Galactic nucleus is 

A/„ = 9 X 10« (M/10«Me)^ ( J|j) ' (^) ' x M, (8) 

where a is the halo core radius and all other variables are as defined above. 

This places a strong limit on the presence and survival of clumps. However, there 
are severe caveats the this model. Once 2 clumps have reached the center, they will 
form a binary, and should eventually coalesce into a single body due to energy loss. 
If however, a third dark matter clump arrives before coalescence occurs, then the 
'slingshot' mechanism can eject one of the clumps and the remaining clumps could 
also escape due to the recoil. Most clumps could escape eventual destruction by this 
method, and hence this mehtod of destruction should not be siginificant. In addition, 
dynamical friction will deplete the number of stars in the nuclues, thereby suppressing 
dynamical friction. 

Thus the cores of the clumps are expected to be stable, although the outer layers 
may stripped off by the effects given above. We conclude that a significant fraction 
of dark matter probably survives in the form of clumps. 

2.4 Collisions 

Lake [|l^] analysed the stability of these dark matter clumps against collisions. The 



timescale for disruption owing to mutual collisions of uniform- density spheres yelds 
the critical value for survival versus destruction as : 

Pc > -"y^ (9) 

^ 'Peril 

where / is the fraction of dark matter in clumps with cluster masses M^ < M x 
10^ Mq, Pc is the density of WIMPs inside the clump, p^ is the density of halo dark 
matter, and Pcrtt is the critical density of the Universe. 

If the density of the clump is larger than this critical value, then collisions are not a 
major mechanism for disruption. This is since the lower density clumps cannot survive 
collisions and are thus destroyed, while the heavier clumps survive such collisions. 
Instead, when the density of the clumps satisfies this inequality, disruption by tidal 
shocking becomes the dominant destruction mechanism of the clumps. However, for 
reasonable parameters, the clumps survive disruption by this mechanism also. Hence, 
clumps are expected to be stable against tidal disruption. 



2.5 Experimental Searches for Dark Matter 

Clumping would have a significant impact on searches for dark matter. The nonde- 
tection of dark matter candidates may be due to the lower value of the local halo 
density. An observer in another part of the galaxy where the local halo density was 
higher would observe a multitude of events that would not be observable at other 
parts |]n| . This holds for direct detection methods. In such methods the dark matter 



is detected by direct interactions of the dark matter particle itself with the particles 
in the detector, analogous to the attempt to 'capturing' it, or observing its trail of 
escape. 

Indirect detection techniques involve the search for signatures of the existence of 
dark matter rather than attempting to directly capture it, such as searching for the 
annihilation products of dark matter particles. Such methods would be facilitated 
by clumping, as a larger number of such products would be observed on Earth in 



the direction of a clump |]14[. For example, the study of the high energy gamma-ray 
particles produced due to the annihilation of dark matter particles will be much more 
frequent in clumps of dark matter, and hence the excess prodution in the direction 
of clumps can be calculated. These are constrained by the results from experiments 
such as EGRET (to be explained below), which can thus place limits on the amount 
of clumping. 

Recently, the neutralino has emerged as a more favourable dark matter candidate. 
The discovery of a possible annual modulation signal in recent experimental searches 



for dark matter are compatible with the neutralino as a dark matter candidate |15 



3 Astrophysical Signatures 

3.1 Dark Clumps 

It has been observed that several collections of stars on the sub-galactic scale are 
dominated by dark matter. This may be taken as already existing experimental 
evidence in favour of clumped dark matter. At present, the results are generally 
approximated by 'clumps' of dark matter within a background uniform halo. Some 
of the recently discovered dark-matter dominated clumps are : 

• Extreme dwarf spheroidal satellites of the Milky Way, Draco and Ursa Minor 
are dominated by dark matter. They are essentially dark matter clumps. Their 
existence is difficult to reconcile with Gaussian initial conditions |T3| . 

• The Sagittarius dwarf spheroidal, discovered in 1994, resides in the Milky Way 
halo, at a distance of about 15 kpc from the centre of the Milky Way. It is 
situated well inside the galactic halo potential and about 23 kpc from the Sun, 



and is dark matter dominated 13 



• R.C.Duncan analyzed the absorption line spectra of the double quasar Q2345+007 
and found it to be lensed by a low-luminosity object at z = 1.49 with a mass 
contained within a volume of lateral extent ~ 40h~^kpc [|TB|. This indicates the 
existence of large clumps of dark matter in interstellar space. 

V. Berezinsky, A. Bottino & G. Mignola [|17| showed that the recent MACHOs 



discovered in the galaxy can be interpreted as dense neutralino objects. However, 



the gamma-ray flux is many orders of magnitude higher than the observed one [|17 
Thus, the onserved gamma-ray flux strongly constrians the fraction of dark matter 
that can have survived in clumps in the case of the neutralino dark matter candidate. 

3.2 Gamma- Ray Sources 

Clumps of dark matter could produce gamma-rays and act as gamma-ray sources 
through the annihilation of WIMPs contained in the clumps of dark matter. Lake also 
showed that the clumps would survive collisions with one another and the disk [J13l. He 



showed how several light sources in the Cos B catalogue of gamma-ray sources could 
be candidates of gamma-ray producing dark matter clumps. He proposed Geminga, 
the second-brightest gamma-ray source above 50 MeV as a clump of dark matter. 

3.3 Gam,m,a-Ray Halo 

There is strong statistical support from EGRET (Energetic Gamma- Ray Experiment 
Telescope) data for the existence of a gamma-ray halo surrounding the galaxy. It has 
been suggested that this halo was the result of annihilations of halo dark matter par- 
ticles. Bergstrom, Edsjo and P.Ullio set forth a model involving moderately clumped 
neutralinos 0, |18|. 

In this model the gamma rays are produced by the pair anihilation of neutralinos 
S , the lightest supersymmetric particle in the Minimal Supersymmetric Standard 
Model (MSSM). They assumed that the average clump mass distribution follows 
the somooth component distribution. It was not possible to reproduce the energy 
spectrum (antiproton flux vs. energy) assuming a smooth dark matter halo. But it 
was compatible with the clumped dark matter secnario. 

A strong excess of gamma rays with energy above 55 MeV was also detected 
towards the galactic centre of the Milky Way, This can be explained by the occurence 
of a higher degree of clumping towards the galactic centre. The degree of clumping 
required implies a measurable excess of antiprotons at low enerigies. There is support 
for this from recent BESS measurements. 

Thus the measured excess of cosmic gamma-rays and antiprotons could be ex- 
plained by invoking the existence of clumped dark matter. Bergstrom et al pointed 
out upcoming experiments would be able to rule out or more strongly confirm this 
theory. 

The integrated 7 - ray flux above an energy threshold Eth is given by 

7 



(t>^ (Eth, An, ^) ~ 1.87 X 10-^5 (Eth) x ( J(^)) (Afi) cm-'^s-^sr-^ (10) 

where S{Eth) is a particle physics dependant term. The dependance of the flux on 
the dark matter distribution is contained in the term {J{'ip)) (Af2). This contribution 
of the clumps to the flux of gamma-rays, a measure ot the contribution of clumps to 
the gamm-ray halo of the galaxy, is given by [l^ : 



^jmAm = v^^,i^f5l,j^ I di [ ^ :;;:,:> (n) 



}_ Lf^ / dn: I dij ^^^'^' 

9>.bkpc/S.Vl J An Jiineof sight y0.3GeV/cm^ 

where /S.VL is the angular acceptance of the detector pointing in a direction which 
forms an angle ip with respect to the galactic center and / is the fraction of dark 
matter concentrated in clumps and p{l-,ip) is the dark matter distribution at the 
point {l,ip)- In the smooth case the dependance is quadratic in the density p. The 
relative strengths of the clumped components depends on both the halo profile and 
the product of the halo fraction in clumps and their overdensity, fS. 

It is observed that there is a discrepancy between the 2.75 power law expected 
from a model of diffuse galactic background and the results obtained by EGRET. 
The dark matter signal carrying an excess of 7 rays from a few GeV to the mass of 
the dark matter particls M^ may explain this discrepancy. In addition the smooth 
dark matter halo models cannot reproduce the EGRET data, but this is possible with 
clumped dark matter models. 

The new high energy cosmic gamma-ray detectors both ground and space-based 



shall have the possibility of searching for dark matter signals |^ . 



3.4 Antiproton Flux 

It has been suggested that the antiproton flux in cosmic rays may be due to annihi- 
lating dark matter. In addition, the antiproton flux is able to place constraints on the 
level of clumping as shown by Mitsui, Maki and Orito |2T[|. The antiproton flux in 



cosmic rays may be due to annihilating neutralino dark matter, evaporating primor- 
dial black holes (PBH's) and superconducting strings. The annihilation of neutralinos 
may produce a detectable flux of antiprotons. This process occurs via the path 

EE —>■ q, leptons and other particles —>■ p + other particles (12) 

Using a standard model in the minimal A^ = 1 supergravity with radiative break- 



ing of the electroweak gauge symmetry. Berezinsky et al p2[ pointed out that the 
expected value of Vt^h"^ lies in the range of 0.2-I-— 0.1 in most cosmological models from 



the viewpoint of the age of the Universe. Thus, Mitsui et al pT| took Vt-^h^ = 0.18. 



Now, neutralinos with this value are most likely pure binos, which predominantly 
annihilate into bottom quark pairs bb or tau lepton pairs r^r^ . The bb pairs can give 
antiprotons, but the tau leptons do not. Utilising a Monte Carlo simulation based 



on the diffusion model utilising source spectra of p ontained from the fragmentation 
functions constructed from JETSET, Mitsui et al found that the annihilation rate 
of neutralinos per unit volume ( 5* ) was proportional to the square of the local halo 
dark matter density p=, but that for PBH's, the relation was linear. 

They also obtained fluxes of antiprotons that are too small if a homogenous, 
isothermal and spherical dark matter distribution of density p = 0.3Gev/cm^. They 
noted three casees, however, where the dark matter density would be enhanced, lead- 
ing to a larger antiproton flux. 

Firstly, the Galactic halo may be flattened towards the galactic disk. This could 
lead to an enhancement of the local halo density by a factor of 2, leading thereby to 
an enhancement in the antiproton flux by a factor of 4. 

Secondly, an NGS (non-dissipative gravitational singularity) may reside at the 
galactic centre, leading thereby to a halo distribution of 

Ph{r) = PhQ{r/rQ)~^-^ for r > O.lpc (13) 

where Ph{r) is the local halo density, r© is the galactic radius of the Sun 
and phQ is the halo density near the Sun. This would lead to an enhanced an- 
tiproton flux, possibly by a factor of 200. 

Thirdly, if the halo consists of clumps of dark matter enhancements of 10^ — 10^ 
may be generated in the early universe. If a few per cent of the neutralino dark matter 
is in such clumps, then the antipoton flux would be enhanced by a factor of 20. 

4 Geophysical Signatures 

Possible geophysical signatures of clumping include such large-scale effects as volcan- 
ism and mass extinctions and the possible creation of defects in mica. 

4.1 Volcanism 

Volcanism is one of the means through which heat generated inside the Earth escapes 
it. The primary sources of heat in the Earth are believed to be radioactivity and 
primordial accretional heat. As per conventional wisdom these are the only major 
sources of heat in the Earth; in the case of special planets like lo other sources like 
tidal heating can become significant. Recently it has been shown by two of the authors 
p3| that dark matter accumulations inside the Earth could produce huge amount of 
heat with drastic consequences. They proposed that the annihilation heat due to 
dark matter could be a new source of heat. 

Here the term 'volcanogenic dark matter' is used to refer to the volcanism arising 
due to the capture and annihilation of dark matter inside the earth. This process 
occurs via several stages: Firstly, the capture of dark matter, followed by the accu- 
mulation of the captured dark matter in the centre of the Earth and the annihilation 



of the same. This then leads to plume formation via the process outlined below. The 
plume, on reaching the surface, leads to massive volcanism of the flood basaltic type. 

4.1.1 Capture of Dark Matter 

The capture of dark matter particles by the Earth and planets was first analysed 



by Press & Spergel ||2^. This capture of dark matter particles was investigated on 
account of the neutrinos and other annihilation products produced. These products 
are currently the subject of ongoing experiments [Q. Krauss, Srednicki and Wilczek 
investigated the neutrinos produced by captured dark matter particles, with scalar or 
Dirac neutrinos as examples |2^. They obtained greatly enhanced capture rates for 



masses greater than 12GeV, and used the luminosity of Uranus to place constraints 
on dark matter candidates. Upgoing muons produced by accumulated dark matter 



inside Earth and Sun have also been the subject of investigation ^7 . 

Subsequently Gould obtained greatly improved formulae with enhanced capture 
rates. The Gould formula for the capture is 



Ne = (4.0 X 10^^ sec-') poA^Q'f U{1 - ^ [,^ )6(^)\ (14) 

where po.4 is the halo WIMP density normalized to OAGeVcm"^ , Q = N — {1 — 
4:si'n?9w)Z ^ N — 0.12AZ, f is the fraction of the Earth's mass due to this el- 
ement, A^ = (3t'^/i)/(2'0^/i_), /i = mx/mN, /i+ = (/U + l)/2, /i_ = (/i — l)/2, 
C,i{A) is a correction factor, v = escape velocity at the shell of Earth material , 
V = SkT^/mx = SOOkms'' is the velocity dispersion, and = v^/vesc^ is the dimen- 
sionless gravitational potential. 

This formula neglected the effect of the finite optical depth of the Earth. If this 
effect is taken into account, multiple collisions occur, and this enhances the capture 
by a factor of 5 — 30 x [^. Since the Earth is in the potential well of the Sun, 
the WIMPs actually move faster, making them more difficult to capture directly. 
Gravitational diffusion however, leads to an increase in the phase-space density of 
free-space WIMPs due to encounters with the Earth, Jupiter and Venus, so that 
direct capture of bound and unbound WIMPs is approximately equal to the case of 
direct capture in free space [^ . This means that the formula given above holds good. 



Gould showed that ndirect capture, namely the capture of WIMPs that are in Solar 



orbits, was insignificant |^. Gould also obtained a high-precision WIMP-capture 
formula for the Sun 



4.1.2 Annihilation of Dark Matter 

This dark matter that is captured gradually drifts towards the centre of the Earth. 
This process occurs because the WIMPs lose energy. Here they accumulate till the 
critical threshold is crossed, beyond which annihilation becomes very important. 
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Above this threshold, virtually all the WIMPs captured are annihilated. In the 
WIMP mass range 15 GeV-100 GeV the Gould formula yields total capture rates 
of the order of lO^^sec"^ to lO^^sec"^ . According to the Gould equation above, this 
yields Qe ~ lO'^lV — lO^^PF for a uniform density distribution. 

Interestingly intriguing signatures of a dark matter WIMP have been recently 
obtained at 56 GeV [^ . This mass of WIMP, nearly the same as that of iron, would 
lead to a resonant enhancement in the capture of WIMPs by Earth. The consequence 
would be an increase in the heat production of the Earth, and this would be in support 
of our model calculations. 

4.1.3 Creation of Plume 

In case of uniformly distributed dark matter, the dark matter heat production is 
much less than the annual geothermic heat output of 4 x lO^^iy. However, during 
the passage of a clump core, the heat production due to annihilating dark matter 
exceeds that due to other sources by several orders of magnitude. What would be the 
consequence of these vast amounts of heat ? 

As per conventional theory, this would lead to the creation of what are referred 
to as 'deep mantle plumes'. As per standard geothermologic models, the lowest layer 
of the mantle ( called the D" layer) absorbs heat from the core, thereby decreasing in 
density. Eventually, once a critical minimum density is crossed, this D" layer breaks 
up into rising plumes. These plumes are an efficient way of heat transfer. Ultimately, 
the largest plumes are capable of reaching the surface. 

4.1.4 Flood Basalt Volcanism and Geomagnetic Reversals 

On arrival at the surface, the plume will melt its way through the crust, leading to 
initial explosive volcanism (due to the molten crust) followed by a much longer period 
of flood basaltic volcanism. As is observed from the vast relics left behind by these 
episodes, such volcanism is known to be the most most extensive form of volcanism 
in the world. Examples are the massive Deccan flood basalt volcanic province, the 
Siberian flood basalt volcanic region, and the Brazilian flood basalt zone. Many 
such flood basalt episodes occur simulatenous to major mass exinctions of life. Much 
evidence of a link between these volcanic episodes and the concordant mass extinctions 
has been established [Q. 

In addition, the large quantities of heat produced inside the core are capable of 
leading to geomagnetic reversals. The rise of temperature inside the core leads to an 
instability in the generation of the magnetic fleld, which may lead to geomagnetic 
reversals. 

A very attractive feature of this model is that it can explain the 30 million year 
periodicity in the record of mass extinctions, the 30 million year periodicity observed 
in the periodicity of flood basalt volcanism, and the 30 million year periodicity cal- 
culated for the crossing of clumps (obtained from the models of clumping). 

11 



4.2 Defects in Mica 

Baltz, Westphal and Snowdon-Ifft studied the effect of dark matter on mica [p2| . 
Tfiey studied two models of tfie density profiles for tfie clumped fraction of the dark 
matter halo, namely the Hernquist profile and the Navarro, French and White (NFW) 
profiles. The Hernquist profile is given by: 



Pir) = 1^^^. (15) 

2tt r[r + a)-^ 

, . GM 

(16) 



r + a 



The NFW profile is given by 



..Ml 
■ «^'-tafl + n (18) 



r \ a 



where a is the scale radius, M the mass parameter, p the density potential, the 
gravitational potential, and r the radius at the point under consideration. 

Now, the rate at which defects accumulate in mica is a function of the angle a 
which the incident particle makes with the cleavage plane. Monte Carlo simulations 
show that the rate at which tracks accumulate, ^ is well approximated by |^ : 



dN(a) , , , , 

— - — = co + ci\cosa\ (19) 

where Cq, Ci are functions of the WIMP mass, the dispersion velocity of the halo 
and the velocity of the earth through the halo. They were calculated assuming the 
mica has remained fixed in position relative to the incident WIMPs. In geophysical 
terms, this corresponds to assuming that the mica has remained geologically station- 
ary, in other words the phenomenon of continental drift has been ignored. They took, 
as an example, a single IO^Mq clump interacting with the earth at different times 
from the formation of the mica, taken to be 440 million years ago, and the present. 
Experiments using neutrons instead of WIMPs as the incident particles conformed to 
this model, and hence it is reasonable. 

Baltz, Westphal and Snowdon-Ifft took for their simulations, a mica age of 440 
million years, which corresponds to two galactic years, ie. two rotations of the earth 
about the galactic centre. For A=10, where A is the mass number, they noticed 
a substantial change in the signal of the clump compared to the background. At 
A=100, the change in the signal is much more substantial. Most of the forward angle 
show a pronounced change in the signal. In the clumped case the signal is much 
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larger and the maximum and minimum directions are very different from the non- 
clumped case. For A=1000, the signal is ten times larger than the non-clumped case 
P^ . The signal from a clump encounter 440 My ago appears to be oriented in the 
same direction for all incident directions. This circumstance arises since the North 
American craton was situated at the equator at this time. For signals occurring 220 
My ago, this degenracy is broken as this continent had by that time moved away 
from the equator. The authors suggested using mica samples from different locations 
from different continents, of different ages and different orientations to obtain a large 
amount of information about the nature of the clumping of the dark matter in the 
halo [0. 
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5 Conclusion 

In addition to arising naturally within cosmological dark matter modles, clumped 
dark matter would be stable on astrophysical time scales and clumps are likely to 
have survived to the present epoch. The existence of clumped dark matter also 
would solve many problems in astrophysics and geology that find either difficult or 
no explanation at all within conventional frameworks. Such cases include the galactic 
gamma-ray halo and the periodicity in terrestrial ffood basalt volcanism. Annihilation 
of clumped dark matter should also be treated as a new source of heat in planetary 
bodies which may have observable consequences. In addition they may have already 
been discovered in the form of dwarf spheroidals. The significance of clumped dark 
matter is such that it calls for further studies. 
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